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Abstract—  The  apparent  constants  (K'|  for  the  ionization  of  carbonic  acid  in  seawater  at  various  salini¬ 
ties  IS,  have  been  fit  to  equations  of  the  form 

In  K;  =  in  Kt  +  AiS'12  +  B,S 

where  K,  is  the  thermodynamic  ionization  constant  in  water.  A,  and  B,  are  adjustable  parameters.  The 
temperature  dependence  (TK)  of  K,,  A,  and  B,  were  of  the  form.  a0  +  a,/T  +  a,  In  T.  Equations  of 
similar  forms  have  been  used  to  analyze  the  ionization  constants  for  water  and  boric  acid  and  the 
solubility  product  of  calcite  in  seawater.  The  effect  of  pressure  on  the  apparent  constants  (KJ'/KP)  have 
been  fit  to  equations  of  the  form  - 

ln<Kf/K“)  =  -(AVP  +  0.5  A/C  P2)/RT  (  J  /  1  >  ^  -  < 

where  the  volume  (AV)  and  compressibility  (AK)  changes  are  polynomial  functions  of  temperature. 
The  equations  generated  for  various  acids  in  seawater  have  been  used  to  examine  the  carbonate  system 
in  seawater.  Equations  relating  the  NBS  and  Tris  pH  scales  have  been  derived  as  well  as  equations 
of  pH  as  a  function  of  temperature  and  pressure.  The  equations  from  Hansson  (1972,  Ph.D.  Thesis, 
University  of  Goteborg.  Sweden)  and  Mehrbach  et  ul.  (1973,  Limnol.  Oceunoyr.  18,  897  907)  have 
been  used  to  examine  the  components  of  the  carbonate  system.  At  a  fixed  total  alkalinity  and  total 
carbon  dioxide,  differences  of  +0.01  m-equiv  kg' 1  in  HC03  and  COj  were  found;  however,  the 
[COj]  and  Pco,  are  nearly  the  same.  The  contribution  of  borate  ion,  B(OH)4  determined  from  the 
equations  of  Hansson  (1972,  Ph.D.  Thesis,  University  of  Goteborg.  Sweden)  and  Lyman  (1957.  Ph.D. 
Thesis.  University  of  California.  Los  Angeles)  differ  by  +0.01  m-equiv  kg' 1  for  waters  with  the  same 
salinity  and  temperature.  •  f7T7^  ,  •'  /  a*  /  •/  '  4.1  <3  r  /VC  V  -/ 
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INTRODUCTION 


Recently  we  had  a  need  to  use  apparent  constants 
for  boric  and  carbonic  acid  in  seawater  to  determine 
the  parameters  of  the  carbonate  system  in  the  Medi¬ 
terranean  Sea  (Millero  et  ttf.,  1978,  1979).  An  exam¬ 
ination  of  the  equations  presently  available  (Edmond 
and  Gieskes,  1970;  Takahaski  et  aL  1970;  Almgren 
et  aL  1975;  Mehrbach  et  al.,  1973)  for  the  salinity 
and  temperature  dependence  of  the  apparent  con¬ 
stants  revealed  a  wide  discrepancy  in  the  form  and 
number  of  parameters  used.  Edmond  and  Gieskes 
( 1970)  have  reviewed  the  earlier  work  of  Lyman  (1957) 
and  Buch  (1932,  1933,  1938,  1951).  In  recent  years 
two  independent  studies  on  the  carbonate  system 
have  appeared  (Hansson,  1972;  Mehrbach  et  al., 
1973).  The  measurements  of  Hansson  (1972)  are 
based  on  a  Tris  buffer  scale  (Hansson,  1973)  while 
the  measurements  of  Mehrbach  et  al.  (1973)  are 
based  on  the  NBS  buffer  scale  (Bates  1964). 

The  purpose  of  the  present  paper  is  to  thoroughly 
analyze  the  measurements  of  these  workers  and  to 
provide  a  statistically  valid  representation  of  the 
apparent  constants  as  a  function  of  salinity,  tempera¬ 
ture  and  pressure. 

EQUATIONS  FOR  APPARENT 
CONSTANTS  AT  ONE 
ATMOSPHERE 

Two  general  criterion  were  used  to  examine  the 
salinity  and  temperature  dependence  of  the  apparent 
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constants  for  the  ionization  of  acids  in  seawater:  (1) 
the  temperature  derivatives  of  the  constants  should 
yield  reliable  thermodynamic  quantities  (AH  and 
ACp,  enthalpies  and  heat  capacities);  (2)  the  equations 
should  be  extrapolated  to  pure  water  when  the 
salinity  goes  to  zero.  As  shown  by  Clark  and  Glew 
(1966)  reliable  thermodynamic  functions  can  be  de¬ 
rived  from  equilibrium  constants  (Kt)  by  using  an 
equation  of  the  form  (TK) 

In  K, ,  =  A  +  BfT  +  C  In  T  +  DT  +  ET2  +  . . .  (1) 

where  A,  B,  C,  etc.,  are  adjustable  parameters.  The 
differentiation  of  eqn  ( 1 )  with  respect  to  temperature 
gives 

AH?  =  RT2(A  in  KJ& T)  = 

R[-B+  CT  +  DV  +  2ET'  +  ...]  (2) 

ACp/  =  (MH,/<5T)  = 

R[C  +  2DT  +  6£T2  +  ...]  (3) 

To  fit  equilibrium  constants  for  the  ionization  of  acids 
over  a  temperature  range  of  50‘  only  the  A.  B  and  C 
terms  are  necessary.  Harned  and  Owen  (1958)  have 
used  a  slightly  different  equation 

In  Kj  =  A'  +  BYT  +  C/T  (4) 

to  fit  the  ionization  constants  of  acids  in  water.  The 
statistical  analysis  of  the  ionization  constants  for  car¬ 
bonic  acid  in  water  shows  that  eqn  (1)  is  slightly 
better  than  eqn  (4)  when  three  parameters  are  used. 
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Table  1.  Parameters  for  the  temperature  dependence  of  the  ionization  of  acids  and  the  solu¬ 
bility  of  CaCOj  in  water* 


In  K't  =  A  +  BfT  +  ClnT 

Acid 

A 

-B 

-C 

a 

Reference 

HjO 

148.9802 

13847.26 

23.6521 

0.0014 

Harned  and  Owen  (1958) 

B(OH)3 

148.0248 

8966.90 

24.4344 

0.0027 

Owen  (1934); 

Manov  et  at.  (1944) 

h2co, 

290.9097 

14554.21 

45.0575 

0.0024 

Harned  and  Davis  (1943) 

HCO] 

207.6548 

1 1843.79 

33.6485 

0.0033 

Harned  and  Scholes  (1941) 

Solid 

Calcite 

303.1308 

13348.09 

48.7537 

0.035 

Jacobson  and  Langmuir  (1974); 
Berner  (1976) 

Aragonite 

303.5363 

13348.09 

48.7537 

0.05 

Berner  (1976) 

*  Temperature  range  0  <  t°C  <  50°. 


We  have  selected  eqn  (1)  to  fit  the  thermodynamic 
ionization  constants  of  water  (Harned  and  Owen, 
1958X  boric  acid  (Owen,  1934:  Manov  et  at.,  1944), 
and  carbonic  acid  (Harned  and  Scholes,  1941; 
Harned  and  Davis,  1943)  in  water  from  0  to  50CC. 
The  values  of  A.  B  and  C  (eqn  1)  along  with  the 
standard  errors  are  given  in  Table  1.  Values  of  AH 
and  ACp  for  the  ionization  of  H20,  boric  and  car¬ 
bonic  acids  at  25  C  determined  from  eqns  (2)  and 
(3)  are  given  in  Table  2  along  with  directly  measured 
values  (Hale  et  at.,  1963:  Vanderzee  and  Swanson, 
1963:  Millero,  1979a;  Larson  and  Hepler,  1969; 
Berg  and  Vanderzee,  1978).  The  agreement  in  the 
measured  and  calculated  values  of  AH'  are  good.  For 
ACp  the  agreement  is  not  as  good  due  to  the  errors 
in  differentiating  eqn  (l)  twice.  As  discussed  by  Lar¬ 
son  and  Hepler  (1969)  a  more  reliable  temperature 
function  for  the  ionization  of  acids  can  be  obtained 
by  integrating  eqn  (3)  than  the  differentiation  of  eqn 
(1).  As  more  reliable  heat  capacity  data  becomes 
available  (for  a  recent  review  see  Millero,  1979a), 
this  should  become  possible  in  the  near  future.  For 
the  analysis  of  seawater  data  our  representation  of 
the  temperature  dependence  of  the  ionization  of  acids 
are  more  than  sufficient. 

We  have  also  fit  the  solubility  product  of  calcite 
(Jacobson  and  Langmuir,  1974  and  Berner,  1976)  in 
water  from  0  to  50  C  to  eqn  (1).  The  coefficients  A,  B 
and  C  and  the  standard  error  are  gi.  m  in  Table  I. 


Berner  (1976)  has  recently  measured  the  solubility  of 
aragonite  in  water  at  5  and  25  C.  His. results  give 
JCSp(aragonite)//Csp(calcite)  =  1.48  ±  0.14  at  25  C  and 
1.43  ±  0.14  at  5°C.  These  values  can  be  compared 
with  the  ratio  calculated  from  the  free  energies  of  the 
two  solids.  Langmuir  (1964)  gives  AGfcalcite  - 
aragonite)  =  -  230  +  30  cal  mol  " 1  at  25  C.  while 
Christ  et  at.  (1974)  give  —218  ±  30 cal  mol  1  at 
25°C  and  —278  ±  40 cal  mol'1  at  5  C.  The  value  of 
K,p(aragonite)/K5p(calcite)  =  1.46  ±  0.08  at  25  C 
obtained  from  these  two  studies  are  in  good  agree¬ 
ment  with  the  measured  value  of  Berner  (1976).  At 
5‘!C  a  value  of  1.65  ±0.12  is  obtained  from  Christ  et 
at.  (1974)  which  is  slightly  higher  than  the  experimen¬ 
tal  value  of  Berner  (1976).  Within  the  error  of  the 
various  estimates  and  measurements,  we  select  the 
value  of  X,p(aragonite)/K,p(calcite)  =  1.5  ±0.1  from 
5  to  25°C.  Combining  this  estimate  with  the  tempera¬ 
ture  coefficients,  we  have  estimated  the  effect  of  tem¬ 
perature  on  the  solubility  of  aragonite  [the  coefficient 
A(aragonite)  =  /((calcite)  +  ln(l.5)  or  0.405]. 

The  equation  selected  to  represent  the  salinity 
dependence  of  apparent  constants  is  given  by 

In  K',  =  A'  +  B'S' 11  +  CS  (5) 

This  equation  was  arrived  at  by  examining  the  value 
of  K\  of  Lyman  (1957),  Hansson  (1972).  Mehrbach  et 
at.  (1973),  Culberson  and  Pytkowicz  (1973),  and  In¬ 
gle  et  at.  (1973,  1975)  at  a  given  temperature.  Figures 


Table  2.  Values  of  AH  and  ACp'  for  the  ionization  of  boric  and  carbonic 
acids  in  water  at  25X 


AH  (kcalmol'1)  ACpTcal mol' 1  K"1) 

This  This 

Acid  study  Literature  study  Literature 


HjO 

13.50 

13.34* 

-47 

-si* 

moH), 

3.34 

3.37' 

-49 

—  721, 

h2co3 

2.23 

2.I9J 

-90 

-72” 

HCOj 

3.60 

3.5  C 

-67 

-51" 

*  Half,  et  at.  (1963);  Vanderzee  and  Swanson  (1963). 
h  Calculated  from  the  heat  capacity  data  tabulated  by  Millero  (1979a). 
‘  Larson  and  Hepler  (1969). 

J  Berg  and  Vanderzee  (1978). 
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Fig.  1.  The  log  of  the  apparent  ionization  constant  of  boric 
acid  (K'b)  in  seawater  plotted  versus  the  square  root  of 
the  salinity. 


1-5  show  the  plots  of  In  X \  vs  S1'2  for  the  ionization 
of  water,  boric  and  carbonic  acids,  and  the  solubility 
of  calcite.  As  is  quite  apparent  from  these  plots  the 
values  of  In  X,  for  boric  acid  and  H2C03  from  the 
work  of  Lyman  and  Mehrbach  are  linear  functions  of 
S'12  and  the  intercepts  approach  the  thermodynamic 
values  for  pure  water.  For  the  ionization  of  H20, 
HCOj ,  the  solubility  of  calcite,  and  all  the  constants 
measured  by  Hansson  (1972)  a  second  degree  func¬ 
tion  is  needed.  The  smooth  curves  in  Figs  1-5  rep¬ 
resent  the  least  squares  best  fit  of  the  data. 

The  temperature  dependence  of  the  terms  in  eqn  (5) 
have  been  determined  by  fitting  the  directly  measured 
apparent  constants  to  the  equation 

In  X;  -  In  K,  =  A'S112  +  FS  (6) 

where  the  parameters  A  and  F  have  the  same  tem¬ 
perature  form  as  In  X,  (the  thermodynamic  constant) 

A  —  flo  -F  t/T  +  In  T  +  . . .  (6a) 

F  —  bo  +  b\fT  +  bi  In  T  +  . . .  (6b) 

The  parameters  a0  to  a,  and  b0  to  bt  for  the  ionization 
of  HjO,  B(OH)j,  H2C03,  HCOj  and  the  solubility 
of  CaC03  are  given  in  Table  3  (along  with  the  stan¬ 
dard  errors).  The  least  squares  computer  program 
used  to  determine  the  adjustable  parameters  was  the 
same  as  used  in  our  equation  of  state  work  (formu¬ 
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lated  by  Manny  Mehr).  Successive  higher  order  terms 
were  dropped  until  the  standard  error  increased  suffi¬ 
ciently  to  indicate  by  an  F-test  that  the  term  was 
significant  at  the  950/0  level.  The  b ,  and  h2  terms  were 
not  needed  for  any  of  the  systems  studied.  The  b0 
term  was  not  needed  for  the  apparent  constants  of 
H2C03  and  B(OH)3  measured  on  the  NBS  pH  scale. 
The  a2  term  was  needed  only  for  the  ionization  of 
water.  It  should  be  pointed  out  that  Hansson's 
measured  values  of  In  K't  and  In  K'w  at  various  tem¬ 
peratures  were  made  only  at  S  =  35“,,,,. 

The  derived  parameters  represent  the  best  statisti¬ 
cally  representation  of  the  experimentally  measured 
apparent  constants.  The  validity  of  the  constants  for 
seawater  diluted  with  pure  water  below  20", „,  needs 
to  be  demonstrated  experimentally. 

The  new  carbonate  results  of  Mehrbach  el  a l. 
(1973)  and  Hansson  (1972)  have  similar  standard 
errors  and  appear  to  be  more  reliable  than  the  earlier 
work  of  Lyman  (1957).  The  results  of  Hansson  (1972) 
for  boric  acid  appear  to  be  more  reliable  than  the 


Fig.  2.  The  log  of  the  apparent  ionization  constant  of  car¬ 
bonic  acid  (X i )  in  seawater  plotted  versus  the  square  root 
of  the  salinity. 
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Fig.  3.  The  log  of  the  apparent  ionization  constant  of  car¬ 
bonic  acid  (K'2)  in  seawater  plotted  versus  the  square  root 
of  the  salinity. 


earlier  work  of  Lyman  (1957).  The  pure  water  results 
of  Hansson  (1972)  and  Culberson  and  Pytkowicz 
(1973)  have  similar  standard  errors. 

The  differences  between  the  values  of  In 
obtained  for  the  ionization  of  acids  by  various 
workers  are  shown  in  Figs  6-8.  The  smoothed  values 
of  In  K\  and  In  K'2  for  carbonic  acid  from  the  results 
of  Lyman  and  Mehrbach  are  within  the  sum  of  the 
standard  errors  of  the  fits  (0.06  for  In  K\  and  0.13 
for  In  K'2)  over  the  entire  temperature  and  salinity 
range  of  the  measurements.  The  differences  in  In  K\ 
at  25°C  are  quite  small  (0.004)  and  are  randomly  dis¬ 
tributed.  For  In  K\  the  diffetences  are  independent 
of  temperature  and  nearly  a  linear  function  of  salinity. 
The  errors  at  high  salinities  are  slightly  larger  (0.21 
at  S  =  40)  than  the  sum  of  the  standard  errors  (0.13). 
From  the  comparisons  with  Hansson’s  work  at  these 
high  salinities,  Lyman's  work  appears  to  be  in  error. 


Fig.  4.  The  log  of  the  apparent  ionization  constant  of  water 
(K'»4  in  seawater  plotted  versus  the  square  root  of  the 
salinity. 


The  differences  between  the  values  of  In  K'w,  In  Kg, 
In  K\  and  In  K\  obtained  by  Hansson  and  other 
workers  at  25  C  are  given  in  Table  4.  The  values  of  In 
K\  and  In  K'2  from  Hansson’s  work  are,  respectively. 
0.32  ±  0.03  and  0.37  +  0.04  higher  than  the  values  of 
Lyman  and  Mehrbach  et  al.  from  5  to  40", salinity. 
This  is  equivalent  to  a  pH  difference  of  0.14-0.16 
which  is  in  good  agreement  with  the  measured  values 
of  0.14-0.16  (Hansson,  1973)  and  is  related  to  the 
differences  in  the  NBS  and  Tris  buffer  scales.  For 
boric  acid  the  differences  in  In  KB  (between  Hansson 
and  Lyman)  are  0.22  ±  0.04  which  is  0.1  lower  than 
expected.  The  recent  measurements  and  calculations 


Fig.  5.  The  log  of  the  apparent  solubility  constant  of  calcite 
f Kip)  in  seawater  plotted  versus  the  square  root  of  the 
salinity. 
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Table  3.  Parameters  for  the  temperature  dependence  of  apparent  constants  in  seawater' 


Solute 

«o 

In  K'  -  In  K, 
a i 

=  ( a0  +  a,/T  +  u2  In  T)S‘  2 
a2  ho10! 

+  b0S 

a 

Reference 

HjO 

-97.9429 

4149.915 

14.8269 

-2.3694 

0.019 

Hansson  (1972) 

-79.2447 

3298.720 

12.0408 

-1.9813 

0.020 

Culberson  and 
Pytkowicz  (1973) 

BfOH), 

0.0473 

49.10 

— - 

— 

0.049 

Lyman  (1957) 

0.5998 

-75.25 

— 

-1.767 

0.008 

Hansson  (1972) 

HjCOj 

-0.1110 

74.02 

- 

— 

0.053 

Lyman  (1957) 

0.0221 

34.02 

— 

-- 

0.007 

Mehrbach  et  al.  (1973) 

0.5709 

-84.25 

— 

-1.632 

0.021 

Hansson  (1972) 

HCOJ 

0.8925 

-89.92 

— 

-1.530 

0.097 

Lyman  (1957) 

0.9805 

-92.65 

-3.294 

0.033 

Mehrbach  et  al.  (1973) 

1.4853 

- 192.69 

-• 

-5.058 

0.042 

Hansson  (1972) 

CaCOj 

1.6233 

- 1 18.64 

-6.999 

0.065 

Ingle  et  al.  (1973) 

Ingle  (1975) 

*  In  the  salinity  range  0  <  S",„  <  40  and  temperature  range  0  <  t  C  <  50  . 


of  Byrne  and  Kester  (1974)  at  25°C  give  a  value  of 
Ain  KB  =  0.27  which  is  slightly  higher  (the  earlier 
work  of  Buch,  1933,  gives  Ain  /C8  =  0. 18).  These 
comparisons  tend  to  indicate  that  the  results  of  Hans- 
son  for  boric  acid  are  too  low  by  0.15  in  In  K'B. 
Further  studies  are  needed  to  clarify  these  differences. 

The  differences  in  In  K'w  of  Hansson  (1972)  and 
Culberson  and  Pytkowicz  (1973)  are  0.04  +  0.01  at 
25°C  from  5  to  40",,,  salinity  which  is  nearly  within 
the  sum  of  the  standard  error  of  the  fits  (0.03). 


Another  method  that  can  be  used  to  examine  the 
reliability  of  the  apparent  constants  for  the  ionization 
of  carbonic  acid  is  to  examine  the  ratio  of  K\/K'2 
(Skirrow,  1975).  Values  of  In  (K\/K'2)  obtained  from 
the  derived  equation  are  shown  at  various  salinities 
and  temperatures  in  Fig.  9.  The  results  of  Lyman  et 
al.  are  in  good  agreement  and  within  the  standard 
error  of  the  measurements.  Larger  differences  in 
in  (ATj/AT^)  are  found  in  the  temperature  dependence 
than  in  the  salinity  dependence. 


Fig.  6.  Comparisons  of  the  apparent  ionization  constants  Fig.  7.  Comparisons  of  the  apparent  ionization  constants 
for  carbonic  acid  obtained  by  Mehrbach  et  al.  (1973)  and  for  carbonic  acid  obtained  by  Hansson  (1972)  and  Miiir- 
Lyman  (1957).  bach  et  al.  (1973). 
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Fig.  8.  Comparisons  of  the  apparent  ionization  constants 
for  boric  acid  and  water  obtained  by  Hansson  (1972). 
Lyman  (1957)  and  Culberson  and  Pytkowicz  (1973). 


Fig.  9.  Comparisons  of  the  ratios  of  K\  and  K'j  for  the 
ionization  of  carbonic  acid  obtained  by  Lyman  (1957). 
Hansson  (1972)  and  Mehrbach  et  al.  (1973). 


The  values  of  K[r  for  the  solubility  of  calcite  in 
seawater  of  Ingle  et  al.  (1973,  1975)  have  been 
selected  because  their  results  cover  a  wide  salinity  and 
temperature  range.  These  results  have  recently  been 
confirmed  at  5  and  25  C  by  the  measurements  of 
Morse  et  al.  (1979).  The  recent  measurements  of 
Morse  et  al.  (1979)  yield  values  for  the  solubility  of 
aragonite  in  seawater  that  are  -  1.5  times  their  value 
for  calcite.  Thus,  the  equations  for  the  effect  of  tem¬ 
perature  and  salinity  on  the  solubility  of  calcite  in 
seawater  can  also  be  used  for  the  solubility  of 
aragonite. 


EQUATIONS  FOR  THE  EFFECT 
OF  PRESSURE  ON  APPARENT 
CONSTANTS 

The  equation  selected  lo  represent  the  effect  of 
pressure  on  the  apparent  ionization  constants  of  acids 
in  seawater  (Culberson  and  Pytkowicz.  1968)  is 

In  (Kf/K?)  =  — (AV,/RT)P  +  (0.5  AK,/RT)P2  (7) 

where  AV,  and  AK,  are  the  volume  and  compressibi¬ 
lity  change  for  the  ionization.  The  values  of  AV,  and 


Table  4.  Comparisons  of  the  values  of  In  K,  obtained  by  various  workers  for  various  acids  in  seawater  at 

25  C 


Ain  Kw 

A  In  K  , 

A  in  K\ 

AlnK2 

S (’’,.,) 

a 

b 

b 

r 

b 

c 

5 

0.049 

0.214 

0.264 

0.259 

0.386 

0.290 

to 

0.059 

0.252 

0.323 

0.318 

0.447 

0.359 

15 

0061 

0.259 

0.347 

0.344 

0.453 

0.391 

20 

0.060 

0,252 

0.353 

0.354 

0.433 

0.404 

25 

0.057 

0.235 

0.349 

0.353 

0.396 

0.405 

30 

0.052 

0.212 

0.336 

0.344 

0.345 

0.398 

35 

0.046 

0.183 

0.318 

0.329 

0  285 

0.384 

40 

0.040 

0.150 

0.295 

0.309 

0.218 

0.365 

0.044  ±  0.007 

0.22  ±  0.04 

0.32  ±  0.03  0.32  ±  0.03 

0.37  ±  0.08  0.37  ±  004 

‘  Hansson  Culberson  and  Pytkowicz. 
b  Hansson  Lyman. 

*  Hansson  Mehrrach  et  al. 
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Table  5.  Coefficients  for  eqn  (7a)  and  (7b)1 


Solute 

«0 

U| 

“ 2 

ajlO3 

bo 

h. 

hz 

o(ln  Xp/X°) 

B(OH)j 

29.48 

-0.295 

-0.1622 

2.608 

2.84 

-0.354 

— 

0.0030 

H,CO, 

25.50 

0.151 

-0.1271 

— 

3.08 

0.578 

-0.0877 

0.0043 

HCOJ 

15.82 

-0.321 

0.0219 

— 

-1.13 

0.314 

0.1475 

0.0042 

CaCCV 

48.76 

— 

-0.5304 

11.76 

— 

-0.3692 

0.054' 

*  Determined  from  the  results  of  Culberson  and  Pytkowicz  (1968). 

b  From  the  work  of  Ingle  (1975).  Valid  only  at  S  =  35  and  for  calcite.  The  value  of  -  AV*  should  be 
2.8  cm/mole  lower  and  AK*  can  be  assumed  to  be  equal  to  AKC. 

'  The  a  at  25C  (0.032)  is  2.3  times  lower  than  the  value  at  2rC  (0.072). 


where  AV,  and  A X,  are  the  volume  and  compressibi¬ 
lity  change  for  the  ionization.  The  values  of  AV,  and 
AX,  have  been  fit  to  equations  of  the  form 

-AV,  =  a0  +  di(S-34.8)  +  a2 1  +  a3i2  (7a) 

-  103  AX,  =  b0  +  bj(S  —  34.8)  +  b2t  +  . . .  (7b) 

Coefficients  of  these  equations  generated  from  the 
results  of  Culberson  and  Pytkowicz  (1968)  are 
given  in  Table  5.  Results  obtained  from  the  work  of 
Ingle  (1975)  for  the  solubility  of  calcite  are  also  given 
in  Table  5. 

As  discussed  elsewhere  (Millero  and  Berner, 
1972;  Ward  and  Millero,  1975)  the  values  of  AV,  for 
the  ionization  of  boric  and  carbonic  acids  determined 
from  molal  volume  measurements  are  in  good  agree¬ 
ment  with  those  generated  from  the  work  of  Culber¬ 
son  and  Pytkowicz  (1968)  (see  Table  6).  Values  of 
AX,  for  the  ionization  of  boric  and  carbonic  acids 
determined  from  the  measurements  of  Culberson 
and  Pytkowicz  (1968)  are  also  in  good  agreement 
(Table  6)  with  the  values  calculated  from  partial  molal 
compressibility  data  (Millero,  1979a).  No  attempt 
was  made  to  force  the  seawater  of  AV,  and  AX,  to  the 
pure  water  values  because  direct  measurements  are 
not  available  at  low  salinities  and  molal  volume  cal¬ 
culations  (Millero  et  ai,  1972;  Ward  and  Millero, 
1974)  indicate  that  AV,  and  AX,  are  strongly  depen¬ 
dent  on  concentration  in  dilute  solutions.  The  coeffi¬ 
cients  a,  and  b,  were  obtained  from  S  =  38.5",,,  and 
13.5  C  results  of  Culberson  and  Pytkowicz  (1968). 

The  values  of  AVC  and  A Kc  for  the  solubility  of 
calcite  in  seawater  are  in  reasonable  agreement  with 
the  values  determined  from  molal  volume  data 
(Millero,  1976,  1979b).  The  value  of 
-AVC  =  31.8  ±  1 .0  cm’  mol  1  determined  for  fora- 
minifera  at  25  C  by  Pytkowicz  and  Fowler  (1967)  is 
-  5.0cm5  mol  1  higher  than  the  value  obtained  by 
Ingle  (1975).  These  differences  are  probably  related 
to  differences  in  the  solid  phases.  The  AV*  for  the 
solubility  of  aragonite  should  be  2.8  cm3  mol " 1  lower 
and  AKa  should  be  approximately  equal  to  AXr 
(CRC  Chemical  Handbook,  1975).  The  value  of  AV* 
determined  in  this  manner  does  not  agree  with  the 
measurements  of  Ingle  (1973)  and  Hawley  and 
Pytkowicz  (1969)  for  oolites  at  2°C.  The  work  of 
Pytkowicz  and  coworkers  yields  a  AV  difference  of 
2.9  cm3  moC  1  for  foraminefera  and  oolites  at  2"C 


while  the  work  of  Ingle  yields  a  difference  of 
10.5  cm3  mol  for  calcite  and  oolites.  Part  of  this  differ¬ 
ence  could  be  due  to  the  larger  temperature  coeffi¬ 
cient  found  by  Ingle  (1975)  for  calcite  compared  with 
the  value  found  by  Hawley  and  Pytkowicz  (1969) 
for  oolites.  Further  measurements  are  needed  to  eluci¬ 
date  these  apparent  discrepancies.  Until  this  is  done 
the  effect  of  pressure  on  the  solubility  of  calcite  and 
aragonite  can  be  estimated  from  the  results  given  in 
Table  5. 


CALCULATION  OF  PARAMETERS 
FOR  THE  CARBONATE  SYSTEM 
IN  SEAWATER 

The  various  components  of  the  carbonate  system  in 
seawater  can  be  characterized  by  measuring  at  least 
two  of  the  following;  pH.  A ZCO,  and  Pto,.  The 
equations  relating  the  selected  parameters  to  the  con¬ 
centration  of  HCOj  .COr.  COJ  =  COj  +  H2COj 
are  given  by  Skirrow  (1975).  The  classical  way  of 
determining  the  components  of  the  carbonate  system 
is  to  make  measurements  of  the  pH  and  total 
titration  alkalinity,  A given  by 

A T  =  [HCO,]  +  2[COM  +  [B(OH)*j 

+  [OH  ]-[H4]  +  £B  (8) 

i 

where  [B(OH)i],  [OH'].  [H4]  are.  respectively,  the 
total  concentrations  of  borate  ion,  hydroxide  ion  and 
the  proton,  and  £B~  represents  other  bases  in  sea- 

i 

water  that  can  acoept  a  proton  (H2PO*  .  PO«  .  etc.. 
Edmond,  1970).  At  the  pH  of  natural  seawater  the 
carbonate  alkalinity  (A, )  can  be  determined  from 

-4c  «  -4,  -  [B(OH)4  ]  =  [HCOj  ]  +  2[COj  ]  (9) 

The  concentration  of  the  borate  ion  can  be  deter¬ 
mined  from 

[B(OH)*3  =  Ks[B]-|/(Kg  +  (!„)  (10) 

where  the  total  boron,  [B]T  =  1.212  x  10  JS  (Cm  - 
kin,  1965)  Kb  is  the  apparent  composition  constant 
for  boric  acid,  and  a„  is  the  apparent  activity  of  the 
proton.  From  the  pH  and  Ac  the  parameter  of  the 
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Table  6.  Comparisons  of  AV  and  A K  for  boric  and  carbonic  acids  in  seawater  at 
2S'  C  determined  from  direct  measurements  and  molal  volume  data 


By. 


Dlstrlhutier/ 


Aval  1  ab  i  \1 

|  Avail  r.iui/'ov 
special 


Solute 

-  AV  (cm1 
Meas 

mol " 1 ) 

Calc 

-  AK 103  (cm3  mol 
Meas 

‘atm  ') 
Calc' 

BtOH), 

27.1 

25.7- 

2.8 

2.9 

h2co3 

22.3 

2t.7b 

0.9 

1.7 

hco3 

15.4 

15.9‘ 

2.4 

3.6 

CaCOj(S) 

35.5 

39.7d 

2.5 

8.2 

•  From  Ward  and  Millf.ro  (1975 1.  A  higher  value  of  -  AV  =  26.4  cmJ  mol  1 
is  obtained  for  0.725  m  NaCl. 

h  From  Millero  and  Berner  (1972)  using  P(H2Ol  =  IX. 1  cm1  mo!  '  instead 
of  the  value  of  18.8  used  in  the  original  paper. 

'From  Millero  and  Berner  (1972)  using  P(CO,  I  =  19.2  cm'  mol  1 
(Duedall.  1972). 

d  From  K  data  of  Millero  (1979a). 


:  From  Millero  (1976). 


carbonate  system  are  determined  from 

[HCOJ  ]  =  Ac/[  1  +  2  Ki/a„] 

(ID 

[COM  =  AcK2l(aH  +  2K2) 

(12) 

[COJ]  =  (Aca„IK\)l(\  +  2  Ki/a„) 

(13) 

ICO,  =  [HCOJ]  +  [COj]  +  [C02] 

=  Affl  +  Oh/K,  +  Aj/ajJAl  +  2K2/uH) 

(14) 


(12).  Since  many  workers  do  not  have  computer  sys¬ 
tems  available,  we  have  generated  values  of  pH  from 
eqn  (15)  as  a  function  of  temperature,  pressure  and 
salinity.  The  values  have  been  generated  using  the 
constants  of  Hansson  and  Mehrbach  et  al.  The  effect 
of  temperature  and  pressure  on  the  pHSbs  are  shown 
in  Fig.  9.  The  effect  of  salinity  on  the  temperature 
coefficients  from  S  ~  30  to  40,,1>,  is  quite  small  and 
can  be  neglected.  The  least  squares  fit  of  the  generated 
values  of  pH  on  the  NBS  and  Tris  scales  are  given  by 


where  the  values  of  Kj  are  the  apparent  constants 
for  carbonic  acid.  The  values  of  K \  and  uH  used  in 
these  equations  are  for  the  in  situ  temperature  and 
pressure  of  sample.  The  effect  of  temperature,  pressure 
and  salinity  on  the  apparent  constants  were  given 
earlier.  The  effect  of  temperature,  pressure  and 
salinity  on  the  pH  needed  to  calculate  in  situ  values 
of  uH  can  be  determined  in  two  ways:  (1)  an  iterative 
computer  technique,  and  (2)  by  using  explicit  function 
of  pH  as  a  function  of  S,  T  and  P.  The  iterative 
method  (Edmond  and  Gieskes,  1970)  requires  the 
solving  of  the  cubic  equation  (Skirrow,  1975). 

«h  +  M  -  H  +  K’giA  -  BtyA 

+  «h [K\K‘dA  -  B  -  1)  +  K\K’2{A  -  2 )]/A 

+  K\K\K'AA  -  B  ~  2M A  =  0  (15) 

which  is  derived  by  combining  eqns  (10)  to  (14).  The 
values  of  A  and  B  in  eqn  (15)  are  given  by 


NSB  Scale 


pH, 

=  pH25  +  4(t  -  25)  +  B(t  -  25)2 

(18) 

I03/1 

=  -9.702  -  2.378 (pH25  -  8) 

+  3.885  (pH25  -  8)2 

(18a) 

10*B 

=  1.123  -  0.003  (pH25  -  8) 

+  0.933 (pH,,  -  8)2 

(18b) 

(a  =  0.002) 

Tris  Scale 

pH,  = 

=  pH2J  +  4(1  -  25)  +  B(t  -  25)' 

(19) 

103A 

=  -9.296  +  32.505(pH2,  -  8) 

+  63.806  (pH25  -  8)2 

(19a) 

10  *B 

=  3.916  +  23.000 (pH2,  -  8) 

+  41.637  (pH25  -  8)2 

(19b) 

(o  =  0.003) 

A  =  /1,/ICOj 

(16) 

B  =  [B]t/XC02. 

(17) 

Values  of  A  and  B  are  first  estimated  from  eqns  (8)  to 
(12)  using  the  measured  aH  and  the  in  situ  apparent 
constants.  Equation  (15)  is  solved  for  uH  using  the 
solutions  for  a  cubic  equation  (CRC  Chemical 
Handbook,  1975)  or  by  iterative  solutions  (Ben- 
Yaakov,  1970;  Almcren  et  at..  197S).  Once  the  in  situ 
aH  is  determined  the  various  parameters  of  the  car¬ 
bonate  system  can  be  determined  from  eqns  (8)  to 


where  pH25  is  the  measured  pH  at  t  =  25  C.  These 
equations  are  valid  from  t  =  0  40  C.  S  =  30  40",,,, 
and  pH2,  =  7.6-8. 2.  The  generated  values  of  pH  have 
also  been  used  to  relate  the  two  scales  at  various 
temperatures 

pHnbs  =  pHx,i»  +  A  +  Bl  (20) 

10M  =  5.93  -  3.6(S  -  35)  (20a) 

I03fl  =  3.381  +  0.058| S  -  35)  (20b) 

The  differences  between  the  two  buffer  scales  are 
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Table  7.  Comparisons  of  the  pHNBS  and  pHTris  scales  in 
seawater 


Temp 

S  =  30 

ApH 

S  =  35 

ApH* 

5  =  40 
ApH 

OC 

0.073 

0.057 

0.038 

5 

0.092 

0.077 

0.059 

10 

0.109 

0.095 

0.079 

15 

0.125 

0.113 

0.099 

20 

0.141 

0.130 

0.117 

25 

0.156 

0.147  (0.159)+ 

0.136 

30 

0.170 

0.162 

0.153 

35 

0.184 

0.178 

0.169 

40 

0.198 

0.192 

0.185 

*  ApH  —  pHNBS  —  pHTns. 
t  Experimentally  measured  by  Hansson. 


given  in  Table  7.  The  calculated  value  of 
pHNBs  -  pH-rm  =  0.147  at  25°C  and  S  =  35  is  in 
reasonable  agreement  with  the  directly  measured 
values  of  0.159  by  Hansson  (1973).  Further  measure¬ 
ments  as  a  function  of  salinity  and  temperature  are 
needed  to  confirm  eqn  (20).  The  effect  of  applied  pres¬ 
sure  on  the  pH  (NBS  scale)  of  seawater  is  a  linear 
function 

pHf  =  pH?  +  AP  (21) 

The  slopes  are  functions  of  t  and  S.  A  least  squares  fit 
of  the  generated  values  of  pH’’  —  pH0  give 

— 103/4  =  0.424  -  0.0048(S  -  35)  -  0.00282t 

—  0.0816  (pH0  —  8)  (21a) 
which  is  valid  from  S  =  32-38  and  t  =  0-25°C. 

A  comparison  of  the  various  parameters  of  the 
carbonate  system  calculated  from  the  apparent  con¬ 
stants  based  on  the  two  pH  scales  are  given  in 
Table  8.  These  calculations  were  made  at  t  =  25°C  and 
S  =  35°m  for  waters  of  the  same  AT  =  2.400  m-equiv 
kg' 1  and  SC02  =  2.111  mmol  kg  1 .  The  differences 
of  K„  determined  from  the  results  of  Lyman  (1957)  and 


Hansson  (1973)  cause  a  difference  in  [B(OH)i]  of 
0.01 1  m-equiv  kg'  ‘.  This  difference  is  larger  than  what 
one  would  expect  (0.004)  from  the  standard  errors  in 
In  Kg.  Until  these  discrepancies  are  clarified,  it  is  not 
possible  to  determine  the  carbonate  alkalinity.  Ac,  to 
better  than  +0.01  m-eqiv  kg'1.  The  differences  in 
[HCOf]  and  [C03]  are  0.010 mmol  kg'1.  For 
waters  with  the  same  Ac  (the  values  in  parentheses  in 
Table  8)  the  difference  in  [HCOJ]  increases  to 
0.019  mmol  kg’  while  the  difference  in  [CO^_]  is 
nearly  the  same  (0.009  mmol  kg  '1).  These  differences 
in  [HCOJ]  and  [C03]  are  about  what  one  would 
expect  from  the  standard  errors  obtained  from  the 
constants  of  Hansson  (1972)  and  Mehrbach  et  al. 
(1973).  The  calculated  values  of  [COJ]  and  Pco, 
using  the  two  sets  of  apparent  constants  are  in  good 
agreement  and  within  the  expected  errors  due  to  the 
standard  errors  of  K\  and  K2  These  calculations 
indicate  that  the  presently  available  apparent  con¬ 
stants  expressed  on  the  NBS  and  Tris  pH  scales 
cause  maximum  differences  of  SAC  =  0.011  m-equiv 
kg'1.  <5[HCOj]  =  0.019  mmol  kg'1.  <5[COi']  = 
0.010  mmol  kg'1,  <5[COJ]  =  0.0001  mmol  kg'  *,  and 
(5PCOj  =  4  x  10'6atm.  It  is  interesting  to  compare 
the  errors  due  to  the  apparent  constants  to  those 
obtained  from  the  experimental  errors  in  pH.  Ar  and 
ZC02.  The  precision  of  the  measurements  of  these 
parameters  are  +0.003  in  pH  (Pytkowicz  el  al.,  1966; 
ZKINO,  1975);  ±0.002  m-equiv  kg’1  in  A1  (Hansson. 
1972;  Almgren  and  Fonselius,  1976)  and  +0  003 
mmol  kg'1  in  IC02  (Hansson,  1972;  Almgren  and 
Fonselius.  1976).  The  errors  of  various  parameters 
of  the  carbonate  system  due  to  experimental  measure¬ 
ments  of  pH  and  AT  or  Ar  and  £C02  are  given  in 
Table  9.  The  uncertainty  in  the  various  carbonate 
parameters  due  to  experimental  errors  are  within  the 
errors  due  to  the  apparent  constants,  except  for  the 
errors  in  PCO;  (due  to  errors  in  XC02).  If  one  set  of 
constants  are  used  the  errors  become  comparable 
except  for  HC03'  and  COj".  To  improve  our  ability 


Table  8.  Comparisons  of  the  parameters  of  the  carbonate  system  determined  at  25  C  and  S  =  35 

using  various  pH  scales* 


Parameter 

NBS  Scale 

Tris  Scale 

A” 

At,  m-equiv  kg' 1 

2.400 

2.400 

0 

ZCOj,  106  mmol  kg’1 

2.111 

2.111 

0 

pH 

8.151 

8.007 

0.144 

B(OH); 

0.095  ±  0.003 

0.084  +  0.001 

0.01 1  +  0.003 

Ac.  m-equiv  kg' 1 

2.305  ±  0.003 

2.316  +  0.001 

0.01 1  ±  0.003 

(2.305  +  0.001) 

[HCOj],  mmol  kg' 1 

1.890  ±  0.011 

1.880  ±  0.014 

0.010  +  0.018 

(1.871  ±  0.018) 

<0.019  +  0.018) 

[CO!'],  mmol  kg  ' 

0.208  ±  0.005 

0.218  ±  0.007 

0.010  +  0.009 

(0.217  +  0.007) 

(0.009  +  0.009) 

[CO'2],  mmol  kg ' 1 

0.01340  ±  0.0002 

0.01347  +  0.0004 

0.0000.3  ±  0.0004 

(0.0130  +  0.0004) 

(0.0001  +  0.0004) 

PCOj,  10‘atm 

472  ±  T 

471  ±  14 

1  ±  16 

(468  ±  14) 

(4  ±  16) 

*  The  errors  are  only  related  to  the  errors  in  the  apparent  constants. 
"The  ±  deviations  have  been  determined  from  (A|  +  A2)I,J. 

1  Calculated  using  the  Henry's  law  a  from  Weiss  (1974). 
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Table  9.  Comparisons  of  the  experimental  errors  in  the  values  of  parameters  for  the  carbonate  system  with  errors  due  to 

the  apparent  constants 


From  pH  and  AT 

From  At 

and  IC02 

Parameter 

ApH* 

A4Tb 

A  ATb 

AI  C02‘ 

From  constants'1 

B(OH  )* ,  m-equiv  kg  “ 1 

A,-,  m-equiv  kg" 1 

0.0005 

— 

0.003 

0.001 

0.011 

0.0005 

0.002 

0.001 

0.001 

0.011 

[HCOj  ],  mmol  kg"1 
[COj "],  mmol  kg" 1 

0.002 

0.002 

0.001 

0.005 

0.010 

0.002 

0.001 

0.001 

0.002 

0.011 

[CO‘2],  mmol  kg" 1 

0.0001 

0.0001 

0.0001 

0.0002 

0.0002 

Pco . .  10b  atm 

4 

1 

3 

7 

4 

ICO, ,  106  mmol  kg  " 1 

0.001 

0.002 

— 

— 

— 

ApH 

— 

— 

0.003 

0.005 

— 

*  Error  due  to  dpH  =  +0.003  (pH  =  8.15  and  4T  =  2.4). 
b  Error  due  to  <5AT  =  +0.002  (pH  =  8.15  and  Ar  =  2.4). 

''  Error  due  to  SZC02  =  +0.003  ( AT  =  2.4  and  IC02  =  2.1 1 1). 

11  The  maximum  errors  related  to  the  uncertainty  of  the  apparent  constants  on  NBS  and  Tris  scales. 


to  calculate  reliable  values  of  HCO3  and  CO3 ",  it  will 
be  necessary  to  make  new  measurements  on  K'2.  The 
experimental  errors  in  determining  IC02  result  in 
larger  errors  in  HCOJ  and  fcoj  than  one  might 
desire.  This  results  in  an  uncertainty  of  -0.005  in  pH 
which  is  slightly  larger  than  the  precision  of  pH 
measurements  (Pytkowicz  et  al.,  1966;Zirino,  1975). 
Until  the  precision  in  the  determination  of  IC02 
obtained  by  a  computer  titration  improves,  pH 
measurements  should  be  made.  This  can  be  done  by 
periodic  calibration  of  the  electrode  system  and  using 
the  initial  emf  reading  in  the  computer  titration 
(Almgren  and  Fonselius,  1976).  Although  the  experi¬ 
mental  precision  of  determining  pH,  AT  and  IC02 
using  computer  assisted  titrations  is  quite  high,  the 
accuracy  is  difficult  to  evaluate.  The  results  of  the 
GEOSECS  intercalibration  work  (Takahashi  er  ai. 
1970)  point  out  the  need  of  standards  in  carbonate 
studies  of  seawater. 

The  saturation  state,  ft,  of  calcite  and  aragonite  in 
seawater  can  be  determined  from 

O  =  [Ca][COj]/K;p  (22) 

where  [Ca]  =  2.938  x  10~4  S  mol  kg" 1  (Millero, 
1979c)  and  is  the  solubility  product  of  calcite 
or  aragonite  at  the  in  situ  t,  P  and  S.  The  values  of 
fit  in  this  study  are  based  on  the  carbonate 
constants  of  Mehrbach  et  al.  (1973).  An  error  of 
+0.005  mmol  kg" 1  in  [C03]2"  yields  an  error  of 
±0.1 1  in  ft,  (2.3%)  and  ±0.08  in  ft,,  (2.5%).  Errors  in 
at  1  atm  yield  errors  of  ±0.32  in  ftc  (6.7%)  and 
±0.21  in  ft,,  (6.6%).  From  these  calculations  it  is 
apparent  that  the  largest  source  of  errors  in  determin¬ 
ing  ft  arc  due  to  uncertainties  in  the  apparent  con¬ 
stants  for  the  solubility  of  CaC03.  The  calculation  of 
ft  in  deep  waters  has  the  added  uncertainty  due  to  the 
effect  of  pressure  on  K,p  The  pressure  coefficients  of 
Ingle  (1975)  cause  an  uncertainty  in  ftc  of  -  ±  1.0% 
per  km.  The  differences  between  the  pressure  coeffi¬ 
cients  of  Pytkowicz  and  Fowler  (1967)  and  Ingle 
(1975)  for  calcite  cause  a  larger  uncertainty  in  ft,  of 
-2.4%  per  km.  Further  measurements  are  needed  to 
clear  up  these  apparent  discrepancies. 


To  summarize,  I  have  attempted  to  provide  a  stat¬ 
istically  valid  representation  of  the  apparent  ioniz¬ 
ation  constants  of  water,  boric  and  carbonic  acids, 
and  the  solubility  of  CaC03  in  seawater  as  a  function 
of  t,  S  and  P.  At  present,  I  feel  that  the  apparent 
constants  on  the  NBS  pH  scale  (Lyman,  1957;  Mehr¬ 
bach  et  al.,  1973;  Ingle  et  al.,  1973,  1975;  Culberson 
and  Pytkowicz,  1968)  should  be  used  to  characterize 
the  parameters  of  the  carbonate  system  in  seawater. 
Further  measurements  for  all  of  the  constants  are 
needed  below  S  =  20%,  as  well  as  measurements  for 
K'b  from  S  =  20  40%,  on  the  NBS  scale.  High  pres¬ 
sure  measurements  on  the  solubility  of  various  forms 
of  CaC03  are  also  needed. 
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